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We analyze the generation of soliton-like solutions in a single-component Fermi gas of neutral 
atoms at zero and finite temperatures with the phase imprinting method. By using both the nu- 
merical and analytical calculations, we find the conditions when the quasisolitons, which apparently 
resamble the properties of solitons in non-linear integrable equations, do exist in a non-interacting 
Fermi gas. We present the results for both spatially homogeneous and trapped cases, and emphasize 
the importance of the Fermi statistics and the absence of the interaction for the existence of such 
solutions. 
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I. INTRODUCTION 

Recent pioneering work by DeMarco and Jin Q on 
realization of quantum degeneracy in a trapped gas of 
fermionic ^''K atoms has triggered the interest in prop- 
erties of cold fermions. As opposed to bosons, fermions 
cannot occupy the same quantum state, i.e., they obey 
the Pauli exclusion principle. Since at low temperatures 
only the s-wave scattering gets important, the spin- 
polarized fermions stop to collide at lower temperatures 
and rethermalization process necessary for evaporation 
breaks down. This difficulty has been overcome by trap- 
ping and cooling the fermions in two hyperfine states, en- 
abling in this way the evaporation by collisions between 
atoms in a different spin states |^ . Other groups have 
achieved quantum degeneracy in a gas of ®Li atoms by 
sympathetic cooling of a mixture of bosonic (^Li) and 
fermionic (^Li) atoms reaching the temperatures down 
to {).2Tp |,|4). 

Despite of experimental difficulties in achieving the 
strong quantum degeneracy in fermionic gases, some the- 
oretical work on degenerate Fermi gases in the normal 
phase at zero and finite temperatures has already been 
done. The static properties of one- and two-component 
system were analysed in Refs. |^ and 0. Collective 
excitations in degenerate Fermi gases in the normal phase 
were discussed within the hydrodynamic approximation 
in Rcf. Q and on the basis of the sum rules in Ref. 
1^. In the absence of s-wave collisions other forces, like 
for instance dipole-dipole interaction, start to play role. 
The fermionic dipoles have been investigated recently, in- 
cluding the analysis of stability conditions (lO|. Another 
aspect of dynamic behavior of a degenerate Fermi gas is 
studied in Ref. where the question about the possi- 
bility of generation of solitons and vortices in a normal 
state of non-interacting Fermi gas is asked. It is shown 
by solving the many-body Schrodinger equation that the 
phase imprinting method is capable of generating the 
solitons and vortices in such a system. An alternative 
method of generating such structures is discussed in Ref. 



So far, the phase imprinting method has been success- 
fully applied to generate solitons in Bose-Einstein con- 
densates ||l^, Although not realized experimentally 
yet, it offers also a rather rare opportunity of control 
of the generation of vortices [l^l- The method con- 
sists in passing an off-resonant laser pulse through the 
appropriately designed absorption plate followed by the 
impinging it on the atomic gas. The laser pulse is short 
and strong enough so that the atomic motion is negli- 
gible during the pulse and the basic effect of light is a 
generation of strong optical potential Vopt (r, t) due to 
the ac-Stark effect. This potential acts for a short time 
Tdun and essentially imprints a desired phase of atomic 
wave functions (^(r) = J^"^"'' Vopt{r,t)dt/h. After this, 
the density and the phase try to adopt to each other, 
which results in generation of structures like solitons and 
vortices, depending on the pattern written in the absorp- 
tion plate. 

In this paper we explore further the possibility of 
generating the quasisolitons in a gas of non-interacting 
Fermi atoms. According to the common knowledge re- 
garding solitons, the nonlinearity of the problem plays es- 
sential role, since it cancels the dispersion and allows for 
a propagation of shape preserving pulses. So, at the first 
glance it is surprising that solitons can exist in a system 
of non-interacting particles described by the many-body 
Schrodinger equation, i.e., the linear equation. However, 
we show that, due to the Fermi statistics, there exists the 
regime of parameters where soliton-like solutions can ex- 
ist in the non-interacting Fermi gas. 

The paper is organized as follows. In Sec. II we solve 
the many-body Schrodinger equation for a system of 
non-interacting fermions undergoing the phase imprint- 
ing, leading to the generation of soliton-like structures. 
Both zero and finite temperatures are considered. Sec. 
Ill offers an analytical approach to the problem of prop- 



agation and decay of solitons, based on the Wigner func- 
tion formalism, for a degenerate as well as a Boltzmann 
gas. In Sec. IV another approach based on analysis 
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of density matrices, is followed in order to discuss the 
Thomas-Fermi approximation and to understand the role 
of dimensionality. Finally, we conclude in Sec. 0. 



II. SOLUTION OF THE MANY-BODY 
SCHRODINGER EQUATION 

A. Zero temperature 

At zero temperature the many-body wave function of 
a system of non-interacting particles (fermions) is given 
by the Slater determinant with the lowest available one- 
particle orbitals occupied. 
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The one-particle orbitals ipN^r)) are orthogo- 

nal and they undergo unitary evolution. Hence, their or- 
thogonality is sustained, both during the phase imprint- 
ing and afterwords. Therefore, the one-particle density 
matrix is at any time given by the following formula 
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piir',r",t) = ^*(^"''*)' (2) 

i=l 



and its diagonal part is just the particle density. 

Evolution of each single-particle orbital is obtained 
from the one-particle Schrodinger equation and can be 
split into two stages. First, the phase imprinting tech- 
nique is used to disturb the system in a way of writing 
a desired phase on it. Assuming fast enough phase im- 
printing, the resulting wave function has the form (this 
procedure can be also performed numerically within the 
finite duration of the imprinting process): 

^k{r) — > f'kir) = fkir) exp{i(j>{x)) , 

for k = l,...,iV, where (f>{x) is the phase imprinted on 
each atom. Throughout this paper we deal with the gen- 
eration of quasisolitons, we assume that the imprinted 
phase depends only on a single coordinate, say 'a;', and 
utilize (j){x) ~ (1 -I- tanh(a;/(iphase))- We consider a Fermi 
gas in a three-dimensional box or harmonic trap char- 
acterized by the frequencies lu^, tOy, and ujz- Both are 
separable in Cartesian coordinates. The one-particle or- 



bitals </3fc(r) are taken as a product ipn} {x) ^puy {y) '^nl {z) 
of eigenvectors of the Hamiltonians of one-dimensional 
boxes or harmonic oscillators. 

In the second stage the system undergoes the free evo- 
lution in three-dimensional space. However, because the 
Hamiltonian separates in coordinates 'a;', 'y', and 'z', the 
time propagation of each orbital can be easily reduced to 
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FIG. 1: Images of the fermionic density (normalized to the 
number of atoms) of N=500 atoms after writing a single phase 
step of 1.0 TT and the width dphaso = 2.5 Af onto a one- 
dimensional non-interacting Fermi gas confined in a box at 
different times: (a) 24.7, (b) 49.3, (c) 98.7, and (d) 148.0 
in Fermi units (h/ep)- Each graph consists of two curves: 
one obtained based on the many-body Schrodinger equation 
(solid line) and the second coming from the one-dimensional 
Thomas-Fermi model (dashed line). 



one-dimensional case: 



where 



-iH: 



and Hx is the Hamiltonian of one-dimensional box or 
harmonic oscillator. The diagonal part of one-particle 
density matrix is given then by the expression: 



(3)/ 



(3) 

In Fig. 1^ we plotted the density profiles of a uniform 
system of 500 atoms after imprinting the phase designed 
in the form 0(x) = 0o(l+tanh(a;/dphase))/2. The charac- 
teristic length (the Fermi length) is given by the formula 
Xf — h/pp — 2/NL, where is the Fermi momentum 
and L is the size of one-dimensional box. For the case of 
Fig. 1^, the Fermi length equals 0.004 L and is 2.5 times 
smaller than the width characterising the jump of the 
imprinted phase. As discussed more precisely in Sec. IH, 
the evolution of the density after phase imprinting de- 
pends strongly on the ratio dphasc/AF between the width 
of the phase step and the Fermi length. 

It is clear from Fig. |l| (as well as Fig. ||) that two qua- 
sisolitons (the bright and the dark one) are generated. 
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TABLE L Speeds of dark and bright quasisolitons (in units of 
speed of sound) generated in a one-dimensional homogeneous 
Fermi gas of 500 atoms after imprinting a single phase step 
of TT and various widths dphaso (in units of the Fermi length) . 
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1.25 


-1.061 


+0.953 


2.50 


-1.044 


-hO.958 


5.00 


-1.033 


-^0.963 


12.5 


-1.014 


+0.981 



They propagate in opposite directions with distinct ve- 
locities. We present in Table || the values of soliton's ve- 
locities for various widths of the phase steps. The speed 
of the bright quasisoliton is always higher than the speed 
of sound whereas dark quasisolitons move slower than 
the sound. However, increasing the ratio dphaso/A_F both 
velocities are getting closer to the speed of sound (p_F/m 
in one-dimensional space), reaching eventually (i.e., in 
the limit dphasc S> Af) the value pf/tti as shown in Sees. 
Ill and IV. In Sec. pv| we also analyze the dphasc 3> Af 



limit by decreasing A p (increasing the number of atoms) 
and solving the equations of the Thomas-Fermi model. 

Fig. ^ shows that going deeper into the regime 
rfphaso S> Af the soliton structures become less pro- 
nounced. This a gain is explained in terms of Wigner 
function in Sec . pTl| and based on Thomas-Fermi ap- 
proach in Sec. [IV| . The condition rfphaso S> Af means 
that the height (or depth) of solitons is getting smaller 
in comparison with unperturbed density. Going in oppo- 
site direction (i.e., dphaso ^ Af) results, as presented in 
Fig. ^, in a change of shape of the structures (dispersion) 
on a time scale of the order of h/ep. 

Another aspect of Fig. |l| is connected with the 
Thomas-Fermi approximation discussed in Sec. |IV. 



Each frame of Fig. ^ shows, in fact, two curves; one of 
them (solid line) is obtained by solving the many-body 
Schrodinger equation, and the second (dashed line) is 
representing the Thomas-Fermi approach. Both curves 
match well which proves the validity of the Thomas- 
Fermi approximation in one-dimensional space for large 
enough number of atoms just as it happens in the case 
of a harmonic trap (see Fig. 1 in Ref. where two- 

and three-dimensional cases are also discussed). 



B. Finite temperatures 

One way of including temperature effects on the soli- 
ton's generation and propagation is to allow for populat- 
ing the single-particle states of higher energy than the 
Fermi energy. We have realized this by considering a 
grand canonical ensemble for a system of non-interacting 
fermions. We generate a number of many-particle con- 
figurations according to the Fermi-Dirac statistics, i.e., 
with the probability of populating a particular one- 
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FIG. 2; Comparison of bright -dark soliton structures for vari- 
ous widths of imprinted phase: 2.5 Af (solid line) and 12.5 Af 
(dashed line). The snapshots are taken at times (a) 24.7, (b) 
49.3, (c) 74.0, and (d) 123.3 in Fermi units. Other parameters 
are the same as in Fig. |l[ 
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FIG. 3: Soliton-like structures in the regime where dphasc 
Af (here, dphasc = 1.25 Af). The snapshots are taken at times 
(a) 24.7 (solid line), (b) 49.3 (dashed line), (c) 98.7 (dot- 
ted line), and (d) 148.0 (dashed-dotted line) in Fermi units. 
Other parameters are the same as in Fig. hi. 



particle state ipk given by exp(/3(efe —/i))-|-l)^-'^, where 
is the state energy, /3 determines the bath temperature 
(via /? = l/fcsT), and is the chemical potential (de- 
pendent on temperature in general). More precisely, for 
each one-particle state (unless the probability of popu- 
lating this state is lower than 10~®) we generate a random 
number from the interval [0, 1]. If the number is less than 
exp(/3(efc — ^))-|-l)~^, the one-particle state is populated, 
otherwise it is left empty. Here, Sk is the energy of the 
state under consideration, the chemical potential 



fi^eF[l-{7Ty3){kBT/eFf]. 



(4) 



and £f equals the Fermi energy corresponding to the 
average number of atoms. The above formula for the 
chemical potential applies for the Fermi gas trapped in 
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a harmonic potential for the temperatures up to 0.55 Tp 

i- 

In such a way, having given the temperature and the 
average number of atoms, the many-particle configura- 
tions are sampled. Then, for each many-particle config- 
uration the phase imprinting is performed followed by a 
free evolution. Finally, the density is calculated as an av- 
erage according to a grand canonical ensemble rules. We 
have investigated in this way the propagation of quasisoli- 
tons for temperatures up to 0.3 Tp, already attainable in 
experiment! 2[ |[ ||. 

In Fig. I we plot evolution of the averaged density 
of the Fermi gas confined in a one-dimensional harmonic 
trap of frequency lo after phase imprinting of a single step 
of 2tt and the width equal to 0.5 osc. units, {^/h/ {mui)). 
Here, the chemical potential is given by the formula 
fi — ep + ksTlogll — exp(— ei?/fcsr)], which is a one- 
dimensional counterpart of formula (0) , the temperature 
equals 0.2 Tp, and the averaging procedure is performed 
over 5000 many-body configurations. Both, bright and 
dark solitons are still present in the system exhibiting ba- 
sically the same properties as at zero temperature. For 
higher temperatures (already 0.3 Tp) the soliton struc- 
tures start to get hidden within the noise originating in 
the averaging procedure (Fig. ||). However, one has to 
remember that rather a particular many-body configura- 
tion is realized in a single experiment, and not an ensem- 
ble of them. Finally, Fig. jo] gives an example that going 
from one- to three-dimensional space, predictions made 
based on lower dimensionality survive. Fig. ^ shows 
density profiles of fermionic gas trapped in a spherically 
symmetric potential (averaged over 1000 configurations 
at 0.1 Tp temperature) and can be compared with corre- 
sponding figure from Ref. ]ll} | at zero temperature since 
the details of phase imprinting are the same. 



III. WIGNER FUNCTION APPROACH 

The Wigner function w{x,p,t) is defined as 

■w{x,p, t) ^ J dy (^^\x + y/2, t)^{x - y/2, tfj exp{~ipy) 

where ?/'^(x,t) and 'tljlx^t) denote atomic creation and 
annihilation operators in the Heisenberg picture, respec- 
tively. It obeys the equation 
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FIG. 4: Evolution of the density distribution of 300 (average 
number) Fermi atoms in one-dimensional harmonic trap at 
temperature 0.2 Tf after imprinting a single phase step of 
2.0 TT and dphaso = 0.5 osc. units. The successive frames 
correspond to moments: (a) 0.1, (b) 0.2, (c) 0.3, and (d) 0.4 
in units of 1/lu. 
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FIG. 5: Comparison of bright -dark soliton structures at dif- 
ferent temperatures: 0.1 Tp (less rugged curve, averaged over 
1000 configurations) and 0.3 Tf (more rugged curve, averaged 
over 10000 configurations). Frames (a) and (b) correspond to 
moments 0.1 and 0.3 in units of 1/lj respectively and other 
parameters are the same as in Fig. ^ 



(dt + —dx - dxVextix)dp) w{x,p,t) = 



(6) 



where Vext (x) is the external harmonic trap potential. It 
immediately follows from the definition, Eq. (^) , that the 
density of particles at time t is 



p{x,t) = J ^^(^^'P'*) 



(7) 



The phase imprinting at t ~ modifies the Wigner 
function in the following way: 

wa{x,p,t = 0) ^ 

wo{x,p) ^ J dy (^4,\x + y/2, t = 0)^l,{x - y/2, t = 0) 
X exp{—ipy/h— i(j){x -f y/2) -f i(j){x — y/2)) , 

(8) 
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FIG. 6: Density profiles of non-interacting fermionic gas con- 
fined in tfiree-dimensional splierically symmetric liarmonic 
trap at different times: (a) O.f, (b) 0.2, (c) 0.3, and (d) 0.4 
in units of l/ui after writing a single phase step of (/>o = l.Ovr 
and (, = 0.2 osc. units. The average number of atoms N=10*, 
the temperature is equal to 0.1 Tf, and 1000 configurations 
were used in the averaging procedure. 



where 4>{x) is the imprinted phase. The function Wq{x,p) 
provides the initial condition for the dynamical evolution 
according to Eq. (||). 

Let us consider first the spatially homogeneous case 
{Vext{x) — 0). Eq. with the initial condition 

w{x,p,t = 0) = wq(x,p) can be easily solved. The solu- 
tion is 



w{x,p,t) ^wa{x~ {p/m)t,p) 



(9) 



and reflects the ballistic motion of particles in the ab- 
sence of an external potential. By using the fact that 
in the spatially homogeneous case the correlation func- 

tion {i^^{x + 2//2, t = 0) ^{x - 2//2, i = 0)) = p^{y), i.e., 

it depends only on '?/', the expression for wq{x,p) can be 
rewritten in the form 



wo{x,p) 



dyPoiv) exp{-ipy/n - i<j>{x + y/2) 
+i(l){x-y/2)) . 



(10) 



The above formulae have a general character, and are 
valid for both degenerate and non-degenerate (Boltz- 
mann) Fermi gas. The difference between these two cases 
is only in the specific form of the function po{y). 

Usually in the experiments, the imprinted phase 4>{x) 
is a function with transition regions of the size of the 
order of dphasc between the regions with constant val- 
ues of the phase. Depending on the ratio between dphasc 
and the characteristic width of the function po{y), one 
has completely different behavior of the solution Eq. (||) 
and, therefore, the density evolution. In the degener- 
ate regime, the characteristic width of po is given by the 
Fermi wavelength Xp, which is inversely proportional to 



the Fermi momentum, \p = h/pp. In the opposite case 
of a Boltzmann gas, the width is given by the thermal 
wavelength At = h/^/2niT. 

We first analyze the case where dphasc ^ \f, At . Un- 
der this condition, one can expand the exponent in Eq. 
(|lO| ) in powers of y, and leave only the lowest, linear in 
y, contribution: 

Wq{x,p) « j dypo{y)exp{-ipy/h- i(j)'{x)y) 

= po(p + 0'(a^)) , (11) 

where po{p) is the Fourier transform of the function po(y) 
and (j)' is the derivative of (j). The density of the particles, 
Eq. (0), can then be written as 



p{x,t) 



dp 

dp 
2Tdi 



Wo{x — {p/ 771)1, p) 

PoiP + H'{x-{p/m)t)) . (12) 



Keeping in mind that 0' < Tr/dphase and, therefore, is 
much smaller than the characteristic momenta, pp and 
Pt ^ \JmT for the degenerate and Boltzmann gas, re- 
spectively, Eq. ( p^ ) can be simplified as: 

p[x,t) ^po + J ^p;,{pW{x-{p/m)t), (13) 

where Po = J (dp / {2TTh))po{p) is the initial constant den- 
sity and Pq is the derivative of po. The function 0'(x) 
is only nonzero in the transition regions of the width 
rfphaso- Therefore, Eq. (|l3| ) has a very clear meaning. It 
represents a superposition of bumps moving with differ- 
ent velocities p/m which are distributed with the weights 
P'oip)- 

For a degenerate Fermi gas at zero temperature one 
has Pq{p) = S{p + pp) — S{p — pp), and, hence, 

1 



p{x,t) po + — [<j)'{x - 
Ztt 



{pp/m)t) - (f>'{x - {pp/m)t)] . 

(14) 

The above formula describes a " soliton-like" motion 
(without broadening and interaction) of bright and dark 
structures in the density, described by the function (p'ix), 
in the opposite directions with the same velocity vp = 
pp/m. The physical explanation of this phenomenon is 
very simple. Under the condition dphase ^ Af, the quan- 
tum state of the system after the phase imprinting con- 
tains excitations, particle-hole pairs, only in the vicinity 
of the Fermi surface. (The energy width of the excita- 
tion distribution is of the order of /i/c?phase ^ Pf-) As 
a result, they have the same velocity, either vp or — fi?, 
and, therefore, the corresponding density profile moves in 
both directions with velocities ±vp, respectively, with- 
out any distortion. To be precise, the broadening in 
this case is of the second order in dphase/^F- To deter- 
mine the characteristic time scale Tpo, one has to keep 
in mind that the width of the distribution of the excita- 
tion momenta after the phase imprinting is of the order 
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of 6p fi/dphasc ^ Pf- Hence, the time tfq can be 
estimated as Tpo ~ {fi/Tp){dpha.sc/^F)^ ^ f^/Tp- 

At temperatures T Tp (but T > TpXp /dphasc, the 
condition when the thermal width of the excitation dis- 
tribution is larger then that after the phase imprinting), 
the situation is very similar, and Eq. ( |l3| ) can be written 
in the form 



p{x,t) 



Po 



ds 
- t{pp 



fp{s,T) [(f>'{x + t{pp 
s)/m)] , 



s)/m)) 
(15) 



where f'p{s,T) = {j)p / AmT) coshT^ {spp /2mT) is re- 
lated to the derivative of the Fermi-Dirac distribution 
and describes broadening of the (5-like peaks at p = 
±Pp in the distribution /Og(p) at finite temperatures. 
As a result, now one has a superposition of bright 
bumps and dark wells in the density profile, moving 
with slightly different velocities distributed with two nar- 
row peaks at ±pp/m with the width of the order of 
mT/pp ~ ppT/Tp <C pp. This results in spreading of 
bumps and wells in the density profile on a time scale 
TpT (?i/T)((iphasc/AF) ~ {Ji/Tp)[Tp /T){dpi,^c/ Xp) » 
h/T, h/Tp (at T = TpXp/dphasc one has Tpx = Tpo). 
Therefore, the time evolution of the gas density af- 
ter a phase imprinting can be described in terms of a 
combination of several quasisolitons (bumps) and quasi- 
darksolitons (wells) which move through each other with- 
out any interaction and, for times t t, decay (broaden) 
very slowly. In Fig. ^ we present the results of nu- 
merical calculations of the time evolution of the density 
and Wigner functions for a spatially homogeneous gas at 
the temperature T = O.OlTi? after the 27r~step phase 
imprinting with dphasc = 160Ai?. The results clearly 
demonstrate (see Fig. 0a) that under these conditions 
the soliton-like structures caused by the phase imprinting 
behave like real solitons, i.e., they propagate without any 
noticeable decay (broadening) on a long time scale. The 
blue regions in Figs. 0b and 0c correspond to negative 
values of the Wigner function, and, therefore, indicate 
the quantum origin of the phenomenon. 

For high temperatures T 3> Tp gas is in the Boltzmann 
regime, and Eq. (03) takes the form 



p{x,t) w Po + 



j ^fsip,T)<l>'{x 



{p/m)t) , (16) 



where fB{p,t) = po y/ 2Tr / mT {p / mT) exp{-p^/2mT) de- 
scribes the distribution with two broad peaks at p = ±pt 
with pt V mT, with the width Sp ^ pr- In this case, 
the decay of the original density profile takes place more 
rapidly than in the case of a degenerate gas, on the time 
scale TB - ih/T){dpk,sc/\T) - Tp^T^ > Tp. The 
corresponding numerical results for the density evolution 
at the temperature T = bTp are shown in Fig. ^. Al- 
though the parameters of the phase imprinting are ex- 
actly the same as for the case of Fig. 0, the temperature 
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FIG. 7: The time evolution of the gas density and the Wigner 
function after the phase imprinting with dphasc = 160 Af at 
temperature T = O.OITf: (a) the density profile at the times 
t = 30iF (the dotted curve) and t — lOOtF (the dashed curve), 
respectively, where tp ~ h/Tp; (b) and (c) the Wigner func- 
tion for the time t = lOOtp, where tp = h/ep. 



(or, in other words, the absence of the sharp Fermi sur- 
face) cause a rapid decay of soliton-like structures. 

The same consideration is applied to the trapped gas 
in the harmonic potential. In this case, due to the pe- 
riodicity of the motion, the essential scale for the time 
evolution is set by the inverse level spacing: t < 2tt/uj, 
where uj is the trap frequency. Therefore, after phase 
imprinting one would observe a quasisoliton behavior 
if 27r/(jj ^ Ti?o, tft, or 2'n/uj <C tb for the de- 
generate or Boltzmann regime, respectively. More ex- 
plicitly, in the degenerate regime one has the follow- 
ing requirement for the temperature T and the width 
^phasc- [fta;/max(T, TFAi;-/dphaso)](dphasc/AF) > 1, or 
rfphasc/^o ^ max(l, T/Ziw-v/iV), where Iq = yjh/mui 
is the oscillator length and N is the number of par- 
ticle. Correspondingly, in the regime of a Boltzmann 
gas, the requirement reads: (?iijj/r)(dphasc/AT) ^ 1, or 

rfphaso/^O > yjTjhuj. 
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FIG. 8: The density of the gas with temperature T = STf 
at time t — 5tF (the dashed curve) and t — SOtp (the dotted 
curve), respectively, after the phase imprinting with dphasc ~ 
160Af. 



In the regime where dphasc ^ Af, At, the excita- 
tions created during the phase imprinting have a very 
broad distribution with the width Sp ^ ^/dphasc ^ 
Pf, Pt in a spatially homogeneous gas, and 6e ~ 
^ih/dphasc)'^ in a trapped one. In this case one would 
expect a very rapid decay of "solitons" on a time scale 
T [fi,/max(rF,T)][(iphasc/niin(AF, At)]^ < h/Tp for 
both spatially homogeneous and trapped Fermi gases. 
The corresponding numerical results for the density and 
Wigner function evolution after the 27r~step phase im- 
printing with dphasc — 2Xp are shown in Figs. 9 (for the 
temperature T = O.OlTp) and 10 (for the temperature 
T — Tp). In the quantum degenerate case (Fig. 9), the 
Wigner function has a well pronounced fringes structure, 
and the corresponding complicated structure in the den- 
sity profile decays rather fast. For the non-degenerate 
gas (Fig. |l^), the Wigner function is smooth everywhere 
except for the narrow region subjected by the phase im- 
printing. In this case, the density profile has no small- 
scale structures, as compared to the degenerate case from 
Fig. 9, but also relaxes very fast to the uniform density 
distribution. 



IV. DENSITY MATRIX FORMALISM 

Another way to discuss the soliton's formation in a 
gas of neutral fermionic atoms is to start from a set of 
equations for reduced density matrices. In particular, the 
equation of motion for the one-particle density matrix 
involves the two-particle density matrix p2 and is given 

by [0 

ih-^^Piin,f2,t) = -A-(v? - v2)pi(fi,f2,i) 

+ / [Vin - r') - V{r2 - r')] ^2(^1, r '; fs, f, t) 




,1=1 1r 



T=0_01T, 




FIG. 9: The time evolution of the gas density and the Wigner 
function after the phase imprinting with dphase = 2Af: (a) the 
gas density at the times t = tp (the dashed curve) t = btp (the 
dotted curve); (b) the Wigner function at the time t — tp; (c) 
the Wigner function at the time t = 5tp. The temperature of 
the gas r = O.OITf. 



+ [Vext in , t) - Vext ir2 , t)] pi (fi , f2 , , 



(17) 



where V{ri — r2) is the two-particle interaction term and 
Vext in t) is the external potential. In the limit ri r2 
Eq. (|l^) leads to the continuity equation 



dpin t) 
dt 



+ V ■W,t)vir,t)] 







(18) 



where the density and velocity fields are defined as fol- 
lows: 

p(f,t) lim pi(fi,f2,t) 

vif,t) = JimjVi-V2)xi(n,r2,<) (19) 

and Xi(^i7^2,0 is the phase of the one-particle density 
matrix. 

One can also rewrite Eq. ([T^ ) introducing the center- 
of-mass iR = (ri -I- n)/^) and the relative position (s = 
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FIG. 10: The same as in Fig. ^ but for the temperature 
T = Tf. 



^1 — ^2) coordinates. By taking the derivative of Eq. 
(|l7) with respect to the coordinate s the hydrodynamic 
Eu er-type equation of motion is obtained in the hmit 
s ^ 



dv{r, t) 
dt 



V-T 



v{r, t) ■ V] v{r, t) 



mp{f, t) 
m p{r, t) m 



(20) 



Here, the kinetic-energy stress tensor T, whose elements 
are given by 



Lki = hm — 

m s-^o oskdsi 



(21) 



depends on ai (ri ,r2,t), the amplitude of one-particle 
density matrix and can be calculated based on a lo- 
cal equilibrium assumption which is the substance of 
Thomas-Fermi approximation. Pint {r, t) , on the other 
hand, represents the force due to fermion-fermion inter- 
action (for spin-polarized atomic fermions the dipole- 



dipole interaction is a good candidate) and is defined as 



The one-particle Wigner function within the Thomas- 
Fermi approximation is given by 

w{f,p) = v{n^[67r'p{r}]y^-p') , 

where rjQ is Heaviside's unit step function and hence the 
one-particle density matrix is calculated as 



Pi{R,s) = 



■w{R,p)e 



ips/ h 



After straightforward algebra one obtains: 



pi(i?,s) 



1 



hpF /Pfs\ , h"^ . (Pfs 



The one-particle density matrix is a real function in this 
case and depends only on the length s = |s| of vector 
s. In the limit s ^ the kinetic-energy stress ten- 
sor T is getting diagonal with (S2/m)(l/30/7r2)(67r2p)5/3 
term at each position. Then, assuming the case of non- 
interacting fermions (spin-polarized atoms at low tem- 
perature), the Eqs. ( |l8|) and (|20| ) become a closed set of 
equations for the density and velocity fields: 



|+V.(pz/).0 

— + u • V w + V 
ot 



= 0, 



(22) 



Note that the Eqs. (H) and (|^) are general and could 
be used also for a bosonic system. For example, for the 
system of bosons described by the many-body wave func- 
tion \I'(fi, rjv) = (p{ri) ■ ... • ^p{rN) (as in the case of 
the Bose-Einstein condensate in the mean-field approx- 
imation) the one-particle density matrix is given by 

Pi{R,s) - ip{R+\^ V*{R-\^ ■ 

Now, the kinetic-energy stress tensor possesses off- 
diagonal elements and leads to the so-called "quantum 
pressure" term in the equation of motion. Assuming con- 
tact interaction between the bosonic atoms dominates, 
the appropriate equations are written as 



dp 
dt 
dv 

m 



V -ipv) =0 



cxt 

m 



(23) 



and are the hydrodynamic ||l^ 
Gross-Pitaevskii equation. 



2™^ y/p 
representation of the 
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It is easy to check that for the spatially homogeneous 
system, the one-dimensional counterpart of Eqs. (122* 
has the following solution 



p{x,t) 
v{x, t) 



pa + 77-[0'(a; + {PF/m)t) - (f>' {x - {pF/m)t)] 
Zn 

-l!—[4>'ix + {pF/m)t) + 4)'{x - {pF/m)t)] 



assuming the slope of the phase being imprinted is much 
less than the unperturbed density {4>' <C po)- This as- 
sumption is equivalent to the condition dphasc 3> A_f 
considered by us in Sec. 



[II and the solution ( |24| ) co- 
incides with expression (14). It also coincides with the 
Thomas-Fermi approximation, valid in one-dimensional 
space when the number of atoms is large enough (Sec. 
II). Again, the solution (p^ ) represents two quasisolitons 
(a bright and a dark one) traveling in opposite directions 
after imprinting (at t = 0) the phase (t>{x). The speed of 
quasisolitons is equal to the speed of sound. A similar 
solution can be found also in two- and three-dimensional 
case 



p{t) — Pa A- a[(f)' [x A- ct) — (j)' {x — ct)] 

(t) = - [0' {x + ct) + (b'ix-ct)], (25) 
2m 

where c equals the speed of sound, a '--^ pj^^ or a pg^'^ 
respectively in two- and three-dimensional space, and 
the phase being imprinted changes only in 'x' direction. 
Again, the solutions ( |25| ) are valid in the limit dphaso ^ 
\f. 

It is clear from Fig. || that in the regime where 
rfphase ^ ^F, the dynamics of fermionic system should be 
discussed rather in terms of soliton-like structures with 
dispersion leading to the fast broadening of pulses. Even 
at zero temperature the phase imprinting excites states 
with energy above Fermi level, depending on the width 
of the phase step dphaso-The deviation of the group ve- 
locity of bright and dark wave packets from the zero- 
temperature sound velocity pF/fn maybe qualitatively 
understood by the following argument. In the spirit of 
the local density approximation one considers local value 
of the Fermi momentum determined by the local value 
of the density. Since the local density is lower at the 
center of dark quasisoliton its velocity is thus reduced. 
Of course, the opposite is true for the bright quasisoli- 
ton. Using Eq. ( |24|) one gets local corrections to pp of 
the form SpF — T^0'(O). This correction agrees with 
numerical results in the order of magnitude. More pre- 
cise estimation of the asymmetry of velocities requires, 
however, taking into account dispersion effects and sys- 
tematic expansion in Af /rfphaso up to higher orders. 

There is an equivalent way of calculating the kinetic- 
energy stress tensor T (^l]), involving the one-particle 
Wigner function 



Tfc; — 



PkPi 



w{r,p)(fp, 



(26) 



which is, however, more suitable for the finite- 
temperature extension of the Thomas-Fermi approach. 
In equilibrium, the Wigner function for a uniform sys- 
tem is just the Fermi-Dirac distribution function /(e) — 
(exp(/3(e— p))-|-l)~^ with e — p^ /2m and /i being the par- 
ticle's kinetic energy and the chemical potential (constant 
for a homogeneous case) respectively. The Thomas- 
Fermi approximation assumes, as usually, that the sys- 
tem is locally uniform and the appropriate equations (in 
one-dimensional form) are written as 



dp 
dt 

dv 

at 



1 d 1 



p^/mdp 



dv 

' dx ' mpdx2TTh J exp(/3(fi - /i)) -f 1 



(27) 



Here, the chemical potential is position-dependent and 
is calculated from the relation 



1 

2Trh 



dp 



exp(/3(|--^)) + l 



(28) 



We consider now the case when the temperature is fi- 
nite but the gas is still strongly degenerate, i.e. T <C 
Tf- By using the Sommerfeld expansion and keeping 
only the first term one can easily find the expression 
p — (27r/i)^/(8m) p^. It turns out that in the regime 
where the imprinted phase changes slowly in comparison 
with unperturbed density, the finite-temperature time- 
dependent Thomas-Fermi model defined by the Eqs. 
(P7|) and (E8) has the following solution 



TO V TO 











t) + (f)'{x 



(29) 



For the temperature T = one has {—df/de) = 
S{e — Ef) and the solution ( p9| ) reduces to the one found 
previously (|lj, ^4|). To verify the ansatz ( ^ ) we put it 



in the set of Eqs. (27) and keep only linear terms in 
(according to dphase S> Xf)- The solvability condition 
has then the form 




(30) 



and is true for low temperatures T <C Tf since then 
the derivative of the Fermi-Dirac distribution function 
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TABLE II: Speeds of dark and bright quasisolitons (in units of 
speed of sound) generated in a one-dimensional homogeneous 
Fermi gas of various number of atoms after imprinting a single 
phase step of tt and the width of 2.5 Xf 



N 


c 


c 


500 


-1.061 


+0.937 


1000 


-1.038 


+0.965 


5000 


-1.010 


+0.993 


8000 


-1.006 


+0.994 



strongly peaks at s = ep, whereas the rest of the inte- 
grand is very small. The solution ( p9[ ) is identical with 

A. Numerical results 

We have solved numerically the set of Eqs. (p2|). To 
this end, we applied the inverse Madelung transformation 
Ip^ to Eqs. ( p^ ) and used the Split-Operator technique. 
Imaginary-time propagation method was chosen to gen- 
erate the ground state of a nonuniform (trapped) Fermi 
gas. Real-time propagation, on the other hand, allows 
to investigate the phase imprinting and developing insta- 
bilities after that. 

First, by imprinting the phase on a one-dimensional 
uniform system we again investigate the regime where 
f^phaso S> Xp, however, this time by increasing the 
number of atoms while keeping the same the width of 
the phase step. Of course, larger number of atoms 
means smaller Fermi length. Increasing the number 
of atoms is easily attainable within Thomas-Fermi ap- 
proach as opposed to the direct solution of the many- 
body Schrodinger equation. Results regarding the speeds 
of bright and dark quasisolitons are presented in Table 
. It is clear that reaching the above mentioned regime 



forces both quasisolitons to move with the same velocity 
equal to the speed of sound. 



V. CONCLUSIONS 

In this paper we have studied in detail the possibility 
of generation of soliton-like structures in ultracold Fermi 
gases using the method of phase imprinting. Generation 
of such structures is feasible with existing experimental 
setups. The soliton-like structures are generated because 
in the process of phase imprinting in the limit in which 
the characteristic length of the imprint dphasc <C Aj? only 
a narrow band of states with momenta close to ±pf are 
excited. The life time of soliton-like structures can be 
quite long even at not extremely low temperatures allow- 
ing for observations. We have presented both numerical 
results and analytic theory explaining the mechanism of 
generation and evolution of soliton-like structures both 
in homogeneous quasi ID cylinders and in quasi ID har- 
monic traps. It will be interesting to generalize the pre- 
sented theory to the case of two-component Fermi gas in 
normal phase. 
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